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We discuss neutrino oscillations in curved spacetime. Our heuristic approach can accomodate mat- 
ter effects and gravitational contributions to neutrino spin precession in the presence of a magnetic 
field. By way of illustration, we perform explicit calculations in the Schwarzschild geometry. In this 
case, gravitational effects on neutrino oscillations are intimately related to the redshift. We discuss 
how spacetime curvature could affect the resonance position and adiabaticity of matter-enhanced 
neutrino flavor conversion. 
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I. INTRODUCTION 

Two decades ago an experimental connection between 
quantum mechanics and gravitation was observed jlj]- 
essentially, a gravitational analog of the Aharanov-Bohm 
effect. These neutron interferometry experiments can 
be well described by including a gravitational poten- 
tial energy in the Hamiltonian of the nonrelativistic 
Schrodinger equation. More recently, the effects of grav- 
itation on another quantum mechanical phenomenon — 
neutrino oscillations — have been discussed by a number 
of authors |-|], some of whose results appear to con- 
flict. In this paper we will provide a simple framework 
for studying neutrino oscillations in curved spacetime. 
We hope to clarify some issues left unclear by previous 
treatments [|]-|4|, and provide a more transparent route 
to some of the results obtained in Ref. § by more formal 
methods than we employ here. 

As an illustration of our treatment, we will do ex- 
plicit calculations in Schwarzschild geometry. For radi- 
ally propagating neutrinos, the oscillation formulae will 
be different from those that appear in flat spacetime. 
This is not surprising, since the gravitational redshift for 
radial neutrino propagation is well known J(|. The ex- 
istence of a circular orbit for massless particles in the 
Schwarzschild geometry allows computation of the oscil- 
lation phase for purely azimuthal propagation as well. In 
this case, it will be seen that the oscillation formulae are 
identical to those obtained in flat spacetime. 

We hope our calculations will clarify some issues that 
we feel have been left unclear by previous treatments. 
The authors of Refs. f|,[§ (who in turn follow Ref. @) use 
a semiclassical approximation that employs the action 
of a massive particle as a quantum phase for each mass 
eigenstate. While such an approach might be appropri- 
ate for counting fringe shifts in interference experiments 
0, its application in Ref. 0] to the calculation of the 
neutrino oscillation phase involves the assumption that 
the mass eigenstates are emitted at different times. The 
physical validity of this approach is suspect. In addition, 
it does not allow for possible gravitational effects on the 
spin of the neutrinos || (see Sec. III). Furthermore, Refs. 



|2|,[7| separate out a "gravitational contribution" to this 
phase. Such a separation is only possible for weak fields, 
making the "gravitationally induced phase" a concept of 
limited utility. Related to this, some authors have used 
methods that appear to mix flat- and curved-spacetime 
thinking . These particular treatments can leave the 
reader confused about the precise meaning of quantities 
such as the "energy," and the nature of the coordinates 
(i.e., Do they reflect proper time and distance?). Covari- 
ant calculations do not suffer from these difficulties of 
interpretation, and are thus preferable. 

In Sec. II we formulate a standard treatment of neu- 
trino oscillations in a more geometric framework. In 
Sec. Ill we generalize our treatment to curved space- 
time, with calculations in Schwarzschild geometry in Sec. 
IV (vacuum oscillations) and Sec. V (Mikheyev-Smirnov- 
Wolfenstein (MSW) effect). Conclusions are given in Sec. 
VI. We set G — h = c = 1 throughout this paper. 

II. SIMPLE GEOMETRIC TREATMENT OF 
NEUTRINO OSCILLATIONS: FLAT SPACETIME 

In this section we briefly review a simple, standard 
treatment of neutrino oscillations j^] , and then present a 
geometric version. 

In a standard treatment, the neutrino state is written 

|* Q (M)> = u «s exp HO®* - p 3 x )] Wi)- (i) 

3 

Here flavor (mass) indices are in greek (latin) letters. 
The matrix elements U a j comprise the transformation 
between the flavor and mass bases. The subscript a on 
the left hand side indicates that the neutrino was in flavor 
state a at the initial position x = and time t = 0. The 
mass eigenstates are taken to be energy eigenstates with 
a common energy the three-momenta of the mass 
eigenstates are then 



For discussion on whether neutrino mass eigenstates should 
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where m, is the rest mass corresponding to mass eigen- 
state \vi). To compute the oscillation probability at po- 
sition x, a massless neutrino trajectory is assumed, i.e. 
x = t: 



where 
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and M is the mass operator. After the mass operator has 
done its work, we have 
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This state is then used to compute the oscillation ampli- 
tude. We note that the assumption of a null trajectory is 
necessary for the observation of oscillations; if the mass 
eigenstates could be measured at different positions (or 
times), the interference pattern would be destroyed. Q 

While the form of Eq. (|l|) may not openly suggest 
it, the neutrino state at a given point in spacetime is 
frame- invariant. The quantities on the right hand side 
of Eq. ([!]) — the transformation coefficients between fla- 
vor and mass bases, the phase, the mass eigenstates — all 
are frame independent quantities. Since the connection 
between quantities in flat and curved spacetime is most 
apparent when expressions are written in a manifestly co- 
variant manner, we introduce a generalized form of Eq. 
(0): 
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In this expression, P is the four-momentum operator that 
generates spacetime translation of the mass eigenstates. 
The quantity p nu \\ = dx/dX is the (null) tangent vector 
to the neutrino's worldline x(X) = [t(X), x(X), y(X),z(X)]; 
X is an affine parameter of the worldline. 

We now show that Eq. ([|) is equivalent to Eq. (jl]), 
by simplifying Eq. (Q) for neutrino propagation in the x 
direction. Let iil denote the argument of the exponential 
in Eq. (@). With P = (E, P x , 0, 0), and employing the 
metric rj^ — diagj— 1, 1, 1, 1], we have 
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where u)j is the phase of the j-th mass eigenstate. Since 
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This phase agrees with that in Eq. (||) (in which xq = 0), 
suggesting that the neutrino state as written in Eq. (||) is 
suitable for calculating the vacuum oscillation amplitude. 
Furthermore, the form of Eq. ^ suggests a straightfor- 
ward generalization to curved spacetime. 

We now review how contributions to effective neutrino 
mass arising from neutrino forward scattering off back- 
ground matter can be included in the above formalism. 
These effects are important because they can give rise 
to, for example, the MSW effect. Our treatment is es- 
sentially that found in Ref. JToj ] . 

As an example, we take neutrino propagation through 
an electron background. In this case the Dirac equation 
can be cast in the form 



[fid^+iAtyPl) + Mf]ll>f =0. 



(10) 



(See Ref. 11 for the convention for the Dirac matrices 7 M 
that we employ.) Here ipf is a column vector of spinors 
of different neutrino flavors, and Mf is the vacuum mass 
matrix in the flavor basis: 



U 







(ii) 



with 



be considered momentum eigenstates, energy eigenstates, or 
neither, see for example Refs. and references therein. 

These technicalities are not crucia in the present context. 

2 One might argue that to account for observation of the in- 
terference pattern, an average (rather than null) trajectory 
should be employed; see Ref. In this case we would have 
x = (P/E)t, where P is the average momentum of the mass 
eigenstates having a common energy E. However, this rea- 
sonable argument has no consequences for our study, since 
it merely introduces an overall phase common to all eigen- 
states. We will here follow (perhaps unfortunate) tradition 
and employ null trajectories throughout this paper. 
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The vector is the flavor-basis effective potential ma- 
trix for interaction with the electron background: 
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In this expression, Gf is the Fermi constant, and N£ — 
n e u^ is the number current of the electron fluid; n e is 
the electron density in the fluid rest frame, and u M is 
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the fluid's 4-velocity. Vl is the left-handed projection 
operator. (An index labeling the neutrino's helicity must 
now be included in the eigenkets of Eq. (0)). The form 
of Eq. (10) suggests that the mass shell relation read 



(P» + Ap> L ){Pp + A f J> L ) = —Mf. 
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This expression can be derived by iteration of the Dirac 
equation, with the assumption that the neutrino momen- 
tum is much larger than the inverse scale height of the 
background matter. Assuming that the electron back- 
ground is at rest with respect to the oscillation experi- 
ment, keeping only terms to first order in Gf, and writing 
P= (E,P X , 0,0), we find 
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The three-momentum operator P x now includes effective 
mass contributions from background matter, and can be 
used in Eq. M. 



III. SIMPLE GEOMETRIC TREATMENT OF 
NEUTRINO OSCILLATIONS: CURVED 
SPACETIME 

In applying Eq. (Q) in curved spacetime, we see that 
evaluation of the argument in the exponential can be- 
come more involved. This complication results from the 
dependence of the metric on position. Additionally, we 
may worry about another gravitational effect: since grav- 
itational fields can cause gyroscopes to precess, perhaps 
gravitational fields can also cause neutrino spin flips B . 

How can effects on spin be incorporated into Eq. (El)? 
Gravitational effects on the spin arise through the "spin 
connection" r7 appearing in the Dirac equation in curved 
spacetime [ fL2[ (we here ignore background matter ef- 
fects): 



[ 7 a e^(^+r /i )+M]^ = 0. 
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In this equation and in the rest of this section, greek in- 
dices refer to general curvilinear coordinates, while the 
latin indices a, 6, c, d refer to locally inertial (Minkowski) 
coordinates. The tetrads connect these sets of coordi- 
nates. The explicit expression for T M is 
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Effects on spin can be incorporated into the three- 
momentum operator (like Eq. (|6|)) in an analogous man- 
ner to background matter effects. 



We must first simplify the Dirac matrix product in the 
spin connection term. It can be shown that 
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where ?y ab is the metric in flat space and e abcd is the (flat 
space) totally antisymmetric tensor, with e 0123 = +1. 
With Eq. (19), the non- vanishing contribution from the 
spin connection is 
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In these equations, (— g) 1 ^ 2 = [dct(g A1I/ )] 1 / 2 , where g^v is 
the metric of curved spacetime. The expression in Eq. 
( ^p| ) treats left- and right-handed states differently. In 
order to group it with terms arising from matter effects, 
we can without physical consequence add a term propor- 
tional to the identity to obtain 
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Proceeding as in the discussion of matter effects in the 
last section, the 3-momentum operator used in neutrino 
oscillation calculations can be computed from the mass 
shell condition 



(^ + A G ^ L )(P» + A%V L ) = -M 2 
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where we have not included background matter effects. 

An important point is that the gravitational contribu- 
tion Aq is proportional to the identity matrix in flavor 
space, and diagonal in spin space. It cannot induce spin 
flips on its own. Therefore, it will not have any observable 
effects unless there are other off-diagonal terms in spin 
space (e.g., from the interaction of a neutrino magnetic 
moment with a magnetic field) 0). 

Another complication in applying Eq. (^) in curved 
spacetime is related to the nature of the neutrino tra- 
jectories. In flat spacetime, the neutrino trajectories 
are straight lines. The propagation can be taken to be 
in one spatial dimension, and the variable of integra- 
tion becomes the spatial variable corresponding to that 
direction of propagation, as in Sec. II. However, the 
neutrino trajectories in curved spacetime are typically 
parametrized curves involving more than one spatial vari- 
able: x(X) = [x°(X),x 1 (X),x 2 {X),x 3 (X)]. 

For general trajectories it therefore may be convenient 
to leave the affinc parameter A as the variable of in- 
tegration, as in Eq. (|J). The tangent vector to the 
null worldline, p = p nu u = dx/dX, can be found from 
the geodesic equation or Hamilton- Jacobi equation. The 
four-momentum operator P can be constructed as fol- 
lows: (1) take the neutrinos to be energy eigenstates, and 
set P° = p°: (2) demand that the three-momenta of P 
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and p be parallel, i.e. P l — p l (\ — e) with i = 1,2, 3; and 
(3) {P^ + A^Vl^ + A^Vl) = -AT 2 , with A" now rep- 
resenting both matter and "spin connection" contribu- 
tions. For relativistic neutrinos (e <C 1), ignoring terms 
of 0(A 2 ) and 0(AM 2 ), and remembering that p is a null 
vector, we find 
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(Here the indices z, j refer to the spacelike general curvi- 
linear coordinates, not locally inertial coordinates.) From 
this it follows that the quantity P ■ p appearing in Eq. 
(Q) is simply 
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It is convenient to define a column vector of flavor am- 
plitudes. For example, for mixing between v e and u T , 

x{X >={(v r \*(\)))- (26) 
Eq. (|j) can be written as a differential equation for x(A), 

(27) 



where the subscript / denotes 'flavor basis.' Eq. ( p7| ) can 
be integrated (numerically if necessary) to yield the neu- 
trino flavor evolution. A similar equation was obtained 
in Ref. by more formal methods. 



IV. NEUTRINO OSCILLATIONS IN 
SCHWARZSCHILD SPACETIME: VACUUM 
OSCILLATIONS 

While Eq. (1271) may be useful for calculating the neu- 
trino flavor evolution for general neutrino trajectories in 
general spacetimes, it does not yield a great deal of physi- 
cal insight. In this and the following section we do exam- 
ple calculations in Schwarzschild geometry. This example 
geometry is simple enough that the oscillation formulae 
can be cast in a form that resembles the flat space case. 

In this section we contrast radially propagating neu- 
trinos with azimuthally propagating neutrinos in order 
to demonstrate how gravity affects vacuum neutrino 
oscillations. The geometry in a spherically symmet- 
ric, static spacetime can be globally represented by the 
Schwarzschild coordinate system {x^} (t,r,6,(f>). We 
can take the Schwarzschild line element, which serves to 
define these coordinates, as 



ds 2 = g^dx^dx" 

= -e 2 ^dt 2 + e 2A ^dr 2 



-r 2 d9 2 + r 2 sin 2 a ^ 2 
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Using the tetrads 
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direct calculation yields Aq — 0. This is perhaps ex- 
pected from spherical symmetry, and is in agreement 
with Ref. ||, in which terms arising from the spin con- 
nection vanish in the Schwarzschild geometry. Since the 
components of the metric are independent of the timelike 
coordinate t, there is a conserved quantity, the timelike 
covariant momentum component P t = — E*. (The rela- 
tion between covariant and contravariant components is 
Pp. = 9nvP u -) We take the neutrino states to be eigen- 
states of this quantity. 

Denoting a differential proper distance at constant t 
by d£, we can write 



dX = d£ ( gij—rr- 
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where we have used the facts that the neutrino trajectory 
is null and that the Schwarzschild metric does not mix 
time and space components. Using this expression for dX 
and Eq. (pa), we obtain in vacuum 
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In this equation, Ei = E Je e~ \ r > is the energy measured 
by a locally inertial observer momentarily at rest in the 
Schwarzschild spacetime (and presumably at rest with 
respect to the "oscillation experiment"). Therefore, the 
integrand in Eq. (31) is formally the same as the corre- 
sponding integrand in flat space. 

For radial propagation, dl = e A ^dr is a differential 
element of proper distance for constant t, 6, cj>, and so 
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We see that unlike the flat space case, the integral in 
terms of proper distance is not trivial, due to the gravi- 
tational redshift of the "local energy" Ei and the radial 
dependence of d£. In this manner, spacetime curvature 
(gravity) makes its impact on the oscillation amplitude. 

Of course, in vacuum above a spherical, static source 
of gravitational mass M. , we have 
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where the Schwarzschild radius is r s = 2M. Then Eq. 
(E52|) is trivially integrated: 
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Again, the coordinate difference (r — ro) does not reflect 
a proper distance. (The proper distance I corresponding 



to the coordinate difference (r 



ro 



)is£ = JI y/g^dr 



Jr e A< - r ' 1 dr.) Likewise, E r does not represent the neu- 
trino energy measured by a locally inertial observer at 
rest at finite radius, but rather the energy of the neu- 
trino measured by such an observer at rest at infinity. 
It is generally not possible to extract a separate "grav- 
itational phase" from this expression; nevertheless, it is 
clear that gravity has an effect on the oscillations of ra- 
dially propagating neutrinos. In the weak field limit one 
could define a "gravitational phase," however. 

In the Schwarzschild spacetime there are circular orbits 
of radius r = R = 3A4 for massless particles. Consider- 
ation of neutrino oscillations in such an orbit (which we 
take to be in the plane defined by sin# = 1) can provide 
insight into gravitational effects in the azimuthal direc- 
tion. 

For azimuthal propagation in this orbit, d£ — Rd(j) 1 
and the local energy Ei — E* (I — r s /R) 1 is constant 
along the neutrino trajectory in contrast to the case of 
radial propagation. From Eq. ( |3f| ) we find 



M 2 

n = -— 



(36) 



This expression involves local energy and proper dis- 
tance, making it precisely the same as the corresponding 
flat space expression. Gravity has no effect here.^j 



V. NEUTRINO OSCILLATIONS IN 
SCHWARZSCHILD SPACETIME: MSW EFFECT 

In this section we study an example of gravitational ef- 
fects on MSW resonant neutrino transformations ||f3].|f4| . 
In particular, in Schwarzschild spacetime we find the res- 
onance position and calculate the adiabaticity parameter 
for a radially propagating two-flavor neutrino system in 
an electron background with monotonically decreasing 
density profile. 

It will here be convenient to write the neutrino evolu- 
tion equation in terms of column vectors of flavor ampli- 
tudes: 



("e|*(r)) 

(VrMr)) 



(37) 



3 We note that a "gravitational phase" for motion transverse 
to the radial direction is given in Ref . j^] , in contrast with our 
result. The discrepancy arises because those authors attempt 
to use a gravitational potential energy instead of employing 
covariant methods. 



For radial propagation we obtain 



X(r)=exp 
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e^M] + V f {r) 



s A < r >dr 



(38) 



In this equation, M 2 is the vacuum mass matrix in the 
flavor basis. The contribution from the background mat- 
ter is 



Vf(r) 



v(r) 




(39) 



UnNg is the locally 



with v(r) — 2\f2GpE*n e and n e 
measured electron density. 

Eq. ( |38| ) can also be written as a Schodinger-like equa- 
tion, 



. dx(r) 
1 dr 



2E^ 



(40) 



where the effective mass matrix in the flavor basis is 

M) = e [AM+*Ml M 2 + e A ^V f (r). (41) 

The mixing angle in matter, 6, is defined in terms of the 
diagonalization of Mi: 



M 2 = U^MfU 
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U = 



— sin f 



sin# 
cos# 



(42) 



(43) 
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The difference of the squares of the neutrino mass eigen- 
values in matter, A = rh 2 — rh 2 , is given in terms of 
A = m\ — m\ as 

A = e A [(w-e' i, Acos26') 2 + (e <i> Asin26 
The mixing angle in matter is given by 

tan 26 = 



(44) 



2*Asin20 



(-v + e*Acos26>)' 



(45) 



The "resonance" occurs for sin 2 29 = 1, where a mass 
level crossing occurs and A is a minimum. The resonance 
condition is v(r) = e*^' A cos 20, or 

^G F n e = - A cos 26 = A. cos 26, (46) 
2e V E* 2Ei 

where Ei = —p- u = i?»e - * is the redshifting "local 
energy" introduced in the last section. Unlike flat space, 
the resonance condition is here determined in part by 
an energy which may be redshifted from the energy the 
neutrino was born with at the production site. 
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After dropping a term proportional to the identity ma- 
trix that yields an overall phase, Eq. ([Io]) can be written 
in the basis of instantaneous mass eigenstates as 



; dx(r) 
dr 



-A/4E* 
idO/dr 



-idO/dr 
A/4£* 



x(r), 



where 
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(48) 



and \v\) and 1 1>2 ) are the instantaneous mass eigenstates. 
The adiabaticity parameter 7(r), which compares the rel- 
ative magnitudes of the diagonal and off-diagonal terms 
in Eq. ((47|), is defined to be 



7(r) 
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AE*\d6/dr\ 
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For 7> 1, the neutrino evolution can be approximated 
by a constant superposition of slowly varying instanta- 
neous mass eigenstates, except for a small probability 
cxp [— : |7(r res )] for one mass eigenstate to jump to the 
other at resonance, where r rGS is the position of the res- 
onance HI. We find 



7( r res) 



A sin" 29 
e-®E* cos26» 



A sin 2 29 
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-A d , / _$ 
e — me n e 
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where A is the matter potential four-vector introduced in 
Sec. I. Thus the adiabaticity of the evolving neutrino (the 
degree to which the "jump probability" is unimportant) 
is affected by the spatial dependence of the metric, which 
appears for example in the determination of the local 
energy. 



VI. CONCLUSION 

We have developed a simple formalism for treating neu- 
trino oscillations in curved spacetime. This formalism 
can accomodate matter effects and gravitational contri- 
butions to neutrino spin precession in the presence of a 
magnetic field. 

We have done explicit calculations in Schwarzschild 
spacetime (without background matter effects) . Our sim- 
ple formalism has verified the result of Ref. § that grav- 
itational contributions to spin precession vanish in spher- 
ically symmetric, static Schwarzschild spacetime. How- 
ever, we have found that the oscillation formulae for ra- 
dially propagating neutrinos are altered by gravity. This 
alteration results from the metrical properties of curved 
spacetime. The basis of the effects we found are closely 
related to the gravitational redshift, in the case of both 



vacuum oscillations and the MSW effect. In contrast, az- 
imuthally propagating neutrinos show no alterations to 
their oscillation formulae in the spherically symmetric, 
static case. 

In applications where strong gravitational effects on 
neutrino oscillations are of possible interest (e.g., super- 
novae), matter effects will generally make vacuum oscilla- 
tions (and therefore gravitational effects on the vacuum 
oscillation phase) unimportant. However, gravitational 
effects on the resonance position and adiabaticity of the 
MSW effect are of potential interest. 

For example, the requirement for successful r-process 
nucleosynthesis in a neutron-rich post-core-bounce su- 
pernova environment has been used to delineate values 
of neutrino mass difference and mixing angle which fa- 
vor/disfavor heavy element production [fTs] , These limits 
arise because the and v T neutrinos emitted from the 
supernova have a higher average energy than the emitted 
v e neutrinos. Therefore, an MSW resonant transforma- 
tion of i>n or v r neutrinos provides a population of higher 
energy v e neutrinos that tend to drive the material out- 
side the nascent neutron star toward less neutron-rich 
conditions. In order to preclude r-process nucleosynthe- 
sis, the MSW transformation must be sufficiently adia- 
batic (conversion efficiency ~ 30%), and must occur be- 
fore the radius where the neutron-to-proton ratio freezes 
out ("weak freeze-out radius"). As we saw in Sec. V, 
gravitation affects both the adiabaticity parameter and 
the position of the resonance. 

Current supernova models indicate that the general 
relativistic effects we consider here are probably not very 
important. However, the equation of state of nuclear 
matter is not well understood, and it may be that dur- 
ing the time frame of interest for nucleosynthesis some 
proto- neutron stars may become very relativistic p6| . If 
this turns out to be the case, we may hazard the follow- 
ing conjectures regarding gravitational effects on limits 
on neutrino mass difference and mixing angle from r- 
process considerations. The MSW transformation may 
become more adiabatic because of the redshifting energy 
appearing in the denominator of the adiabaticity param- 
eter. This would extend the limit on neutrino mixing to 
exclude smaller values of the vacuum mixing angle. On 
the other hand, the smallest neutrino mass difference ex- 
cluded by r-process considerations is, roughly, that mass 
difference for which the resonance position of an aver- 
age energy /i (or t) neutrino coincides with the weak 
freeze-out radius. The redshifting energy appearing in 
the resonance condition will tend to pull the resonance 
position closer to the neutron star. However, the weak 
freczc-out radius is determined by the competition be- 
tween the weak interaction rates and the expansion rate 
of the material outside the the neutron star, and the 
weak interaction rates are proportional to the square of 
the redshifting neutrino energy. Therefore, for a given 
mass difference, the redshift will reduce the separation 
of the resonance position and the weak freeze-out radius. 
This would weaken the mass difference boundary of the 



G 



excluded region. Of course, these conjectures are prelim- 
inary, as they are based principally on redshift effects. 
Other effects of a supernova core of sufficiently small ra- 
dius for gravitational effects to become interesting may 
be relevant. Such effects might include changes in the 
density scale height and expansion rate of the background 
matter in the supernova envelope, and alteration of the 
neutrino spectrum. In fact, if the v e and V e are signifi- 
cantly redshifted, then r-process nucleosynthesis may be 
precluded anyway Q. 
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